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The heavy quark effects in deep-inelastic scattering in the asymptotic regime Q 2 m 2 can be 
described by heavy flavor operator matrix elements. Complete analytic expressions for these 
objects are currently known to NLO. We present first results for fixed moments at NNLO. This 
involves a recalculation of fixed moments of the corresponding NNLO anomalous dimensions, 
which we thereby confirm. 
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1. Introduction 

The double-differential cross section for unpolarized deep inelastic scattering via single photon 
exchange can be expressed in terms of two structure functions F2 j i{x, Q 2 ). Especially in the region 
of smaller values of x, these structure functions contain large cc-contributions of up to 20-40 %, 
denoted by F 2 " L (x,2 2 ) l . The perturbative heavy flavor Wilson coefficients H2^{z,Q 2 ,m 2 ) corre- 
sponding to these structure functions are known at NLO semi-analytically in z-space, |JIJ], with a 
fast numerical implementation in Mellin-A^ space given in Oj. Due to the size of the heavy flavor 
corrections, it is necessary to extend the description of these contributions to 0(a^), a s := CC S / (4jt), 
and thus to the same level which has been reached for the massless Wilson coefficients, [j3|]. This 
will allow for a more precise determination of partem distribution functions and the QCD-scale A. 
A calculation of these quantities in the whole kinematic range at N N LO seems to be out of reach 
at present. However, in the limit of large virtualities Q 2 , Q 2 > 10m 2 in the case of F 2 fC (^, Q 2 ), one 
observes that Fj^i, Q 2 ) are very well described by their asymptotic expressions, neglecting 
power corrections in m 2 /Q 2 . In this kinematic range, one can calculate the heavy flavor Wilson co- 
efficients analytically. This has been done for F£ c (x, Q 2 ) to 2-loop order in [f|,|5|] and for F[ c (x, Q 2 ) 
to 3-loop order in In the latter case, the asymptotic result becomes valid only at much higher 
values of Q 2 . As first steps towards the 3-loop calculation at asymptotic scales for F 2 cc (x,<2 2 ) we 
calculated the 0(£) terms of the 2-loop heavy operator matrix elements (OMEs), [Q, which con- 
tribute to the asymptotic 3-loop heavy flavor Wilson coefficients via renormalization. In the present 
paper, we report on new results concerning moments of the heavy OMEs A^$^ q , A qq q and A gq Q 
at 3-loops. In doing so, the fermionic terms of the even moments N = 2. ..12 of the corresponding 
N N LO anomalous dimensions given in [§] are confirmed in an independent calculation. 

2. Heavy Flavor Operator Matrix Elements 

In the following, we consider massive OMEs of the flavor-decomposed twist-2 operators between 
partonic states 

2 °° 2 

Here, S and NS are the singlet and non-singlet contributions, respectively, i denotes the outer on- 
shell particle (i = q,g) and stands for the quarkonic {k = q) or gluonic {k = g) operator emerging 
in the light-cone expansion. The subscript H indicates that we require the presence of heavy quarks 
of one type with mass m, while pL is the renormalization scale. We work in Mellin-space, with the 
Mellin-variable denoted by N. The logarithmic terms in m j\l 2 are completely determined by 
renormalization and contain contributions of the anomalous dimensions of the twist-2 operators. 
Thus at NNLO the fermionic parts of the 3-loop anomalous dimensions calculated in Refs. ^ 
appear. All pole terms provide a check on our calculation. The single pole term allows for a first 
independent calculation of the terms <=c Tp in the 3-loop anomalous dimensions. 
As outlined in Ref. [Qj, in the limit Q 2 S> m 2 one applies the massive renormalization group equa- 
tion to obtain the factorization of the heavy flavor Wilson coefficients into a Mellin-convolution of 

'We consider extrinsic heavy flavor production only. 
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the light flavor Wilson coefficients with the corresponding heavy OMEs. The light flavor Wilson 
coefficients are known up to three loop order [^] and carry all the process dependence, whereas the 
heavy quark OMEs are universal, process-independent objects and contain all mass corrections to 
Hi,l up to terms proportional to (m 2 /Q 2 ) a , a > 1. 

A related application of the heavy OMEs is given when using a variable flavor number scheme to 
describe partem densities including massive partons. The OMEs are then the transition functions 
when going from n/ to n/ + 1 flavors and thus one may define parton densities for massive quarks, 



see e.g. Ref. [10]. This is of particular interest for heavy quark induced processes at the LHC, such 



as cs — > W + at large enough scales Q 2 . In this context, one may show as well that various sum 



rules follow from momentum conservation, QIOQ, e.g: 

,NS,+ 



A qq.Q 



+4 PS 



N=2 m 



= . (2.2) 

N=2 



Note, that in Eq. (2.2) two different PS contributions can be distinguished and we adopt the notation 



^ P Q{q)q := ^2? +A^ S g, where for the {Qq}- term the operator couples to the heavy quark and for 
the {qq, <2}-term to a light quark. 

3. Renormalization 

We work in Feynman-gauge and use dimensional regularization in D = 4 + e dimensions, applying 
the MS-scheme, if not stated otherwise. Renormalization proceeds in four steps, which we will 
briefly sketch here and refer to [Q, ^ [7| ||] for more details. Mass renormalization is performed 



in the on-shell scheme [11], whereas for charge renormalization we use the MS-scheme. The 
remaining singularities are of the ultraviolet and collinear type. The former are renormalized via 
the operator Z-f actors, whereas the latter are removed via mass factorization through the transition 
functions T. After coupling- and mass renormalization, the renormalized heavy flavor OMEs are 
then obtained via 

A=Z- l AT~ l , (3.1) 
where quantities with a hat are unrenormalized. Note, that in the singlet case operator mixing 



occurs and hence Eq. (3.1) should be read as a matrix equation, contrary to the NS-case. The Z- 
and T-factors can be expressed in terms of the anomalous dimensions of the twist-2 operators in 
all orders in the strong coupling constant a s , cf. [J7|] up to 0(a 3 ). For this purpose, we adopt the 
convention 

y = Hd\nZ{}i)/dn . (3.2) 



From Eqs. ( p.l|j3.2| ) one can then infer that for operator renormalization and mass factorization at 
0(a 3 s ), the anomalous dimensions up to NNLO, together with the 1-loop heavy flavor OMEs 
up to 0(e 2 ) and the 2-loop heavy OMEs up to 0(e) are needed. Higher orders in e enter since they 
multiply Z— and T-f actors containing poles in e. This has been worked out in some detail in Ref. 
[f7|], where we presented the 0(e) terms a^, a^g S and flgj PS in the unpolarized case. The terms 

~®ggQ an( ^ ~®mQ were gi ven i n Refs. Thus all terms needed for the renormalization at 3-loops 
in the unpolarized case are known by now. 
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4. Moments of the Massive OMEs at three Loops 

The heavy flavor OMEs at 0{o? s ) are given by 3-loop self-energy type diagrams, which contain a 
local operator insertion. The external massless particles are on-shell. The heavy quark mass sets 
the scale and the spin of the local operator is given by the Mellin-variable N. The steps for the 



calculation are the following: We use QGRAF, [12], for the generation of diagrams. Approxi 



mately 100 diagrams contribute for each case at hand. For the calculation of the color factors we 



use [13]. After undoing the contraction of the operators with the light-like vector A, (A 2 = 0), the 
diagrams are genuinely given as tensor integrals. Applying a projector then provides the results for 
the diagrams for the specific (even) Mellin moment under consideration. The diagrams are further 



translated into a form, which is suitable for the program MATAD [14], through which the expansion 



in £ is performed and the corresponding massive three-loop tadpole-type diagrams are calculated. 



We have implemented all these steps into a FORM-program, cf. [15], and use TFORM, [16], for 
parts of the calculation. We checked our procedures against various complete two-loop results and 
certain scalar 3-loop integrals and found full agreement. 

5. Results 



Applying Eq. ( |3.1| ), one can predict the pole structure of the unrenormalized results and thus the 
logarithmic terms of the renormalized OMEs. Since we consider only terms involving at least one 
heavy quark, we adopt the following notation for the anomalous dimensions 

f=Y(n f +l)-Y(n f ). (5.1) 

As an example, we show the structure of the renormalized result in the PS case, where all renor- 
malization constants are taken at nj flavors. 

= -^+4(» /+ l)ft, fi + 6ft}ln=(^) + {^((„ /+1 )ft, e -ft) 

-(0) (0) -(1) "1 ? ( 42), PS (0) 

+£((»/ + 1 - ®) - *?- y ($) + { 0- - - ¥«§ 

-(0) (0) , ! "I 2 

r Yqg Ygq (JO) JO),./ , UR ■ fi o.\ , n f+ l JO) (2) I , (m \ (0) (2) 

-M j^— [Ygg - Yqq +4(«/ + l)P0,Q + OPo J H — Jcjg d gqQ M n ) + %9 U Qg 

(0) f(0) y 

+4(/3 + Po, Q )a%r +Cf(-<* + \ C2)i 0) # " 47^' PS + I2a^] 
-(n f + l)tg ] ^l Q + a^ S q . (5.2) 

Here, the terms a\ , denote the constant terms in e of the 2-loop OMEs Ajj and p, are the expansion 
coefficients of the j3-function. The subscript Q refers to contributions due to heavy quarks only, 
cf. [01 15[ ^j], and Q is the Riemann ^-function at values i. All quantities in Eq. (5.2) are known for 



(3) 

general values of N, except for ^qL^ which is the genuine 3-loop contribution and remains to be 
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calculated. Note, that it is not possible to factor out nj ■ + 1, not even in the triple pole term. This 
is due to the interplay of the prescription for coupling constant renormalization we have adopted, 
cf. [[1 0], and the fact that the transition functions T apply to sub-graphs containing massless lines 
only. 

We have calculated the OMEs AjSJ, A^ S e + (3) and for N = 2. ..12, using MATAD. All pole 
terms agree with the general structure derived from renormalization. As an example, we show the 
constant term after renormalization of the PS-OME for N = 2 



N= 

A Q{q)q 



/2048 86480% 2 / 3584 53 144 x 2 / 64 128 2 



1280,, 830 \„„ / 128 „ 128,., 9536,, 95638 \ , 

H £s H )C F C A T F + B 4 H G C3 H Cp7> . (5.3) 

27 w 2187/ V 9 5 81 s 729 / F 

The moments we obtain for N = 10, 12 for the corresponding 3-loop anomalous dimensions are 
shown in Table [[]. These results, as well as the moments we obtain for lower values of N, agree with 
the results of Refs. A check for the constant terms is provided by the sum rule in Eq. (|2~2|), 
which is obeyed. Additionally, we find for all moments that the terms proportional to £2 disappear 
after renormalization, which is a general observation made in many D — > 4 calculations. The term 



B4 in Eq. ( |5.3[ ) is given by 



B 4 = -4C 2 ln 2 2 + ^ln 4 2-yC4 + 16Li 4 (^) , (5.4) 

and appears in all OMEs we calculated. Since it does not appear in the light-flavor Wilson coeffi- 
cients, cf. [J3j] , it occurs as a genuine mass effect. 

Table 1: Results for jV = 10, 12 for the 3-loop anomalous dimensions 

_N J^/Tf/Cf 

1218139408 7168 / \ 18846629176433 529979902254031 

10 T>(l + 2nf) + [C A -C F )( 3 C A + C F 

363862125 1 495 V / 11767301122500 323600780868750 

13454024393417 5056 / \ 64190493078139789 1401404001326440151 

T F (1 + 2n f ) + [C A -C F )( 3 C A + C F 

J 429 V / 48885219979596000 3495293228541114000 
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5222779912350 ' 429 V / 48885219979596000 3495293228541114000 

^^■ + /Tf/cV 

27995901056887 192880 / \ 9007773127403 75522073210471127 



131745667845011220000 



192880 / \ 9007773127403 75522073210471127 

T F (l + 2n f )+ [C F -C A )C 3 C A C F 

1 693 V / 97250422500 307518802668000 

13549568 / \ 25478252190337435009 

T F (l + 2 n/ )+ [C F -C a )(s C A 

' 45045 V / 263228107582440000 



3290351344780500 45045 \ 263228107582440000 

35346062280941906036867 



265847305072 50176 / \ 1028766412107043 839864254987192 

10 T F (l + 2 n/ )+ [C F -C A )C3 C A + C F 

420260754375 J 27225 V / 1294403123475000 485401171303125 

2566080055386457 49928 / \ 69697489543846494691 

T F (1 + 2n, + [Cf - C A k 3 C A 



5703275664286200 ' 39039 V 83039693672007072000 

86033255402443256197 



5 1806197823524667520 
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6. Conclusions and Outlook 

We calculated the heavy OMEs A^ S,+ , A^ 3 ] ,PS and A^ Q for even Mellin-moments N = 2... 12 
using MATAD and showed first results. This confirms for the first time in an independent calculation 
the moments of the fermionic parts of the corresponding 3-loop anomalous dimensions, [Q]. We 
expect results for the remaining terms A^g and Agj in the near future, thus enabling us to calculate 
fixed moments of the heavy flavor Wilson coefficients in the asymptotic limit Q 2 ^>m 2 . 

Acknowledgments. The work was supported in part by SFB-TR-9: Computergestutzte 
Theoretische Teilchenphysik and EU-network HEPTOOLS. We would like to thank 
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